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ABSTRACT 

f~^ ■ The photon frequency distribution that results from single Compton scattering of monochromatic radiation on 

(^ ' thermal electrons is derived in the mildly relativistic limit. Algebraic expressions are given for (1) the photon 

(^ , redistribution function, K{v,^ — > v\^'), and (2) the spectrum produced in the case of isotropic incident radia- 

CN ' tion, P{v — !■ v'). The former is a good approximation for electron temperatures kT^ <^ 25 keV and photon energies 

^ |. hv < 50 keV, and the latter is applicable when hv{hv/m^c^) < kT^ < 25 keV, hv <; 50 keV. Both formulae can be 

D ' used for describing the profiles of X-ray and low-frequency lines upon scattering in hot, optically thin plasmas, 

C^ I such as present in clusters of galaxies, in the coronae of accretion disks in X-ray binaries and active galactic 

neclei (AGNs), during supernova explosions, etc. Both formulae can also be employed as the kernels of the corre- 
sponding integral kinetic equations (direction-dependent and isotropic) in the general problem of Comptonization 
on thermal electrons. The Kiiy, f! -^ i^' ,^') kernel, in particular, is applicable to the problem of induced Comp- 
el ■ ton interaction of anisotropic low-frequency radiation of high brightness temperature with free electrons in the 
^ ' vicinity of powerful radiosources and masers. Fokker-Planck-type expansion (up to fourth order) of the integral 
^D , kinetic equation with the P(i> -^ v') kernel derived here leads to a generalization of the Kompaneets equation. 
We also present (1) a simpler kernel that is necessary and sufficient to derive the Kompaneets equation and (2) 
an expression for the angular function for Compton scattering in a hot plasma, which includes temperature and 
photon energy corrections to the Rayleigh angular function. 

Q^ \ Subject headings: accretion, accretion disks — cosmic microwave background — intergalactic medium — line: 

Q> ■ profiles — masers — radiative transfer 
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1. INTRODUCTION 

A monochromatic spectral line will be broadened after single Compton scattering on thermal electrons in an optically thin, hot 
plasma. The emergent spectrum will depend on the angle between the direction, fl, from which the photons are supplied and the 
observer's direction, fl'. A classical problem then arises of finding the redistribution function, K {v,il -^ v' ,i l'), which gives the 



^ probability for the incident photon {v, fl) to be scattered in the direction Jl' with a frequency of v' (Dirac 1925). 
• • In the case in which the incident radiation is isotropic and monochromatic, the spectrum resulting from single Compton scattering 

_ ^ can be found by integrating Kiv, il -^ v' ,^') over the scattering angle (^^ = nn'): 

H; P(i^-^iy')= K(v,n-^i^\n')dn' = 2T: K(iy,n-^iy',n')dfi,. (1) 

Spectra described by equation (IT]) may in particular form when spherical symmetry is present in the system. Consider, for example, 
an isotropic source of monochromatic radiation that is either located at the center of a spherical cluster of galaxies with hot inter- 
galactic gas or distributed spherically symmetrically across the cluster The frequency distribution of photons that have experienced 
a single scattering within the cluster will then be precisely P{i/ -^ v') if the entire cluster is probed at the same time. If instead the 
spectrum is collected from a part of the cluster, then one must take into account the angular distribution of incident photons and make 
use of the direction-dependent K{v, il -^ v' ,^') function. 

Single-scattering line profiles, which can provide important information on the conditions (particularly the temperature) of the 
plasma in the source, can originate in various astrophysical environments. For the X-ray spectral band, these include intracluster gas, 
the coronae of accretion disks around black holes in binary stellar systems and active galactic nuclei (AGNs), and plasma streams 
outflowing from the neutron star during super- Eddington X-ray bursts in bursters. 

High-quality spectroscopic X-ray observations necessary for the detection and measurement of such lines will soon become 
possible with the satellites Chandra, XMM (both akeady in orbit), and Spectrum-X-Gamma, and in a more distant outlook, with 
Constellation-X and XEUS. This was one of the primary motives for us in initiating the current study, in which we calculate analyti- 
cally the functions Kiv, il — > v' ,il') (eq. [R]) and Piv — > v') (eqs. JI^ and [pl[|) in the mildly relativistic limit. 

1.1. Integral Kinetic Equation 

Apart from the single-scattering problem outlined above, Kiiy, il -^ i'' ,il') and Pif — > i'') can be used as the kernels of the corre- 
sponding integrodifferential kinetic equations describing the Comptonization of photons on Maxwellian electrons. In the anisotropic 

1 



2 SAZONOV & SUNYAEV 

problem, the kinetic equation in the case of an infinite homogeneous medium can be written as 

diiiy,n,T) 



dr 



■=- I(iy,n,T)K(iy,n^ iy',n')[l+n{iy',n',T)]diy'dn' 
- / ^liiy' ,n' ,T)Kiiy' ,n' ^ iy,n)[l+n{iy,n,T)]di^'dn' , 



(2) 



where /(j^, fl) is the specific intensity of the radiation, n = c^I /{2hv^) is the occupation number in photon phase space, and r = a-jNeCt 
is the dimensionless time (ctj is the Thomson scattering cross section and A^e is the number density of electrons). It is easy to write 
down kinetic equations similar to equation (||) for any of the standard problems of radiation transfer 
In the isotropic case, the kinetic equation becomes 



dr 



- I(iy,T)P(iy~^iy')[l+n(iy',T)]diy'+ / ^I(iy' ,T)P(iy' -^ i^)[l+n(i^,T)]diy'. 



(3) 



The integrodifferential equation (^ can in general be solved numerically if its ker nel, P(v ~> ly'), is known. Alternatively, one can 
treat Comptonization problems using Monte Carlo methods (see, e.g., the review by Pozdnyakov et al. 1983 ). In the nonrelativistic 
limit, i.e. when the typical photon energy, hi^, and the plasma temperature, kTg, are both negligibly small compared to the electron 
rest energy, njeC^, the variation in intensity at a given frequency is largely governed by transitions in a narrow interval of a continuum 
spectrum near this frequency. If the initial distribution of the photons in frequency is smooth enough (in a problem with a small 
number of scatterings) or the formation of a spectrum by multiple scatterings is studied, it is possible to perform a Fokker-Planck- 
type expansion of the integral equation (g), thereby reducing it to a much simpler differential equation that describes the diffusion 
and flow of the photons in frequency space: 
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The moments of the kernel that enter this equation are found from the formula 

((Az/)")= fp(v-^v')(v'-v)"dv'. 



(4) 



(5) 



Substituting the first two moments, (Av) and {(Av)^), calculated to an accuracy o f kT^/m^c^ and hv/ msC^ (using the kernel given 



by eq. 023] below) into the corresponding terms of equation (W) leads to the famous Kompaneets (1957) equation: 
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(6) 



The Kompaneets equation is valid in the nonrelativistic limit {hv, kT^ <C m^c^). The last parenthesized term in equation (pt) describes 
the frequency diffusion of photons due to the Doppler effect and the transfer of energy from the electrons to the radiation; the first 
term describes the downward photon flow along the frequency axis due to Compton recoil, and the second term, which also is owing 
to recoil, accounts for induced Compton scattering. 

In studying the interaction between energetic photons and hot electrons {hv^ kT^ ^ Q.QlnigC^), relativistic corrections to the kernel 
and the Kompaneets equation become important. In order to find the main corrections (of the order of hv/m^c^ and kT^/m^c^), it 
proves necessary to take into account in equation (Q) all t erms that depend on th e first four moments of the k ernel. The resultant 
generalization of the Kompaneets equation was found by Itoh et al. (1998) and [Challinor & Lasenby (1998 ) when these authors 
examined the distortion o f the spectrum o f the cosmic microwave background (CMB) during interaction with the hot intergalactic 



gas in clusters of galaxies. |^ephaeli (1995 ) pointed out the important rol e of corrections of the order of fcrp/mpC^ in this phe nomenon 



The kernel P{v -^ v') was not given in explicit form in the derivation of [[toh et al. (1998 ) and Challinor & Lasenby (1998). Instead, 
the Fokker-Planck operator was applied to the kinetic equation written in the most general form, where the amplitude appears for 
a transition from the initial state w ith g iven 4-momenta of th e photon and electron into the final state with the corresponding 4- 
momenta. Earlier, Ross et al. (1978) and Qlarionov et al. (1979) added to the Kompaneets equation a dispersion term associated with 
Compton recoil and a term that takes into account, in a first approximation, the transition from the Thomson cross section to the 
Klein-Nishina one. Their equation, which describes much better than the Kompaneets equation the scatterin g of energetic pho tons 
{hv ^ Q.lm^c^) by sufficien tly cold electrons {kT^ <C hv), is a particular case of the more general formula of [toh et al. (1998) and 
Chaflinor & Lasenby (19981). 
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Despite the attractiveness of employing the Fokker-Planck approximation in treating Comptonization problems, its scope is limited. 
For example, in considering the effect of Comptonization from a small number of scatterings on the profiles of narrow spectral lines, 
the initial radiation spectrum cannot be represented as a finite Taylor series in terms of the frequency variation, and, therefore, Fokker- 
Planck-type equations are not applicable. For this kind of problem, it is necessary to make use of the integral kinetic equations (g) or 
(^, which requires knowledge of their kernels, K{v, fl ^ i^' , f!') and P{v -^ v'). 

1.2. The Kernel 



Dirac (1925) has given an approximate expression for the direction-dependent kernel, K{v,^ — > v' ,il'). The Doppler shift was 
taken into account to within the first order in v/c, but Compton recoil was totally neglected. Integration of Dirac's K{v^ fl -^ v' ,il') 
func tion over the scattering angle results in the zero-order approximation for the isotropic kernel P{v -^ v') (H ummer & Mihalas 
1967;see also Weymann 1970). Being symmetric in frequency variation [P{v -^ v + /S.v) = Piu -^ i^-Aiy)], this kernel does not take 
into account the average Doppler increment in the photon energy: {Aiy/i') = AkT^/m^c^. With the zero-order approximation for the 
kernel, it is not possible to describe many important astrophysical phenomena, such as: 

(a) distortion of the spectrum of the CMB in the direction of galaxy clusters ( [Sunyaev & Zeldovich 1972 ;see the recent review by 
Birkinshaw 1 999), 

(b) The y (Zeldovich & Sunyaev 1969) and Bose-Einstein /i ( [Sunyaev & Zeldovich 1970 ) distortions of the CMB spectrum 
resulting from energy release in the early universe (see the reviews by Danese & de Zotti 1977; Sunyaev & Zeldovich 1980; the 
book by Peebles 1993; and the strict constraints on these distortions placed by the Far-Infrared Absolute Spectrophotometer on board 
COBE, Fixsen et al. 1996), 

(c ) the formation of hard power-law tails in the emission spectra of the famous X-ray source Cygnus X-1 ( Sunyaev & Truemper 
1979) an d other stellar-mass black hole candidates, AGNs, and quasars; and in the spectra of accreting neutron stars (e.g., S hapiro et 
al. 1976; ISunyaev & Titarchuk 198q; pozdnyakov et al. 1983t [Tanaka & Shibazaki 1996t poutanen & Svensson 19961; N arayan 
etal. 1998; [Zdziarski et al. 1998^ 

Babuel-Peyrissac & Rouvillois (1970, hereafter BR70) have obtained a more accurate expression for the K{v^ fl — ^ v' , Jl') kernel, 
by taking into account not only the average Doppler and recoil frequency shifts (the latter is {Av/v) = -hv/m^c^), but also, in the 
first order, Klein-Nishina corrections to the scattering cross section. BR70 also attempted to find relativistic corrections due to high 
electron velocities {kT^ ^ O.lnjeC^), but did not include all relevant terms (as we show in the present paper), and therefore their 
expression is val id for nonrelativistic e lectrons only. The allowance for Compton recoil is also made in the K{v, Jl -^ v\^') kernel 



written down by Zeldovich et al. (1972|). Integration of the BR7 kernel over the scattering angle ignoring recoil { hv ^ kT^ <C WgC^) 



leads to the first- order approximation for the P{v -^ v') kernel (Sunyaev 1980). The formula of Sunyaev (1980), as opposed to the 
approximation of Hummer & Mihalas (1967), accounts for the asymmetry of the scattered line profile and enables us to describe (in 
the nonrelativistic limit) the above-mentioned phenomena related to the transfer of energy from hot electrons to radiation. Using this 
expression, one can derive the diffusion part of the Kompaneets operator (the last term of eq. [pfl). The remaining two (flow) terms 
can be found by applying the thermodynamic principle that the desired equation must obey: a Planckian distribution of photons must 
remain unchanged during the interaction with electrons of the same temperature. This is the method followed by the authors of the 
original derivation of the Kompaneets equation for finding not only the term describing the downward motion of the photons along 
the frequency axis, but also the term describing the effect of induced scattering (report No. 336 of the Chemical Physics Institute of 
USSR Academy of Scie nces 1950; Kompaneets 1957; Ya.B. Zeldovich 1968, private communication). 

Kershaw et al. (ivso) have derived a semianalytical relativistic formula for the K{v,^ -^ i^' ,^') kernel that is valid for arbitrary 
values of photon energy and electron tem perature. With this for mula, the calculation of the kernel is reduced to computing a single 
integral over the electron Lorentz factor. Kershaw et al. (1986) further succeded in expanding their basic expression in powers of 



kTe/nieC , and thus derived an algebraic formula for K(y, ft -^ J^', fl') that is applicable in the low-temperature limit kT^ <J OAmgC 
with no limitation imposed on the photon energy. 



Other previous relativistic studies of the Compton scattering kernel also need to be mentioned. Aharonian & Atoyan (1981) were 
the first to derive the exa ct relativistic formula for the direction-dependent kernel, /r(j/,fJ —> ly^Q,'), for an isotropic ensemble of 
monoenergetic electrons. Nagirner & Poutanen (1994) (see also references therein) derived the exact relativistic expression for the 
isotropic kernel, P{i/ -^ v), again for an ensemble of monoenergetic electrons. These formulae make it possible to solve efficiently 
very general Comptonization problems, i.e., with no constraints set on the parameters. However, their application always implies 
the need to perform one or two numerical integrations, usually over a specified distribution of electron energies (e.g., Maxwellian), 
and sometimes over the scattering angle. We also note that the scattering kernel has been the subject of a number of s tudies that 
employed numeri cal methods (e.g., Pomraning 1973|;|lllarionov et al. 197^; pozdnyakov et al. 1979t |Loeb et al. 199l];M olnar & 
Birkinshaw 1999). 

In the present paper, we have obtained algebraic approximate expressions for the kernels Kiv^Q, -^ v' ,^') (eq. |M|) and P(v -^ 
v') (eqs. [|l9|] and [pi]]), which take into account (1) relativistic effects due to high electron velocities and (2) quantum effects, 
namely, Compton recoil and Klein-Nishina corrections. The formula for the angle-dependent kernel, K(i^,n — > v' ,0,'), is a very 
good approximation in the parameter range kT^ <^ 25 keV and hi/ <^ 50 keV. Moreover, in the case of nonrelativistic electrons, this 
formula without the temperature correction terms (eq. |]lO|]) gives the exact result for arbitrary photon energies, including /ij/ 3> mgC^. 
We derived our expression in the same manner as BR70 derived theirs, but our expression is more accurate, since it includes the 
relativistic corrections consistently. 

Our expression for Kii/, il ^ i^' , Q,') can be compared with the corresponding formula derived by Kershaw et al. (1986 ) (their eq. 



[4 1]). Both formulae are fully algebraic and valid in roughly the same electron temperature range {kT^/m^c^ <^ 0.1). The expression 
of Kershaw et al. (198q ) is more general than our equation (0). It holds for any hu in all its temperature terms, while in the case of 
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our formula, this stat ement is true for the m ain, nonrelativistic temperature term only. Therefore, our expression should be derivable 
from the formula of Kershaw et al. (1986). The referee of the present paper informed us that both formulae indeed yield similar 



numerical results for the same parameter values. On the other hand, our expression is significantly simpler in structure and, more 
importantly, when K{v, fl —> v\^') is written in this form, it can be integrated analytically (under some additional constraints on the 
parameters) over the scattering angle or the emergent photon energy. The first of these integrations leads to the algebraic formula 
for the P{v -^ v') kernel, while the second leads to an expression for the angular scattering function. We stress that it is this unique 
integrability that makes equation (0) of interest and useful for further applications. 

We derived the expression for P(v -^ v') by integrating Kiy.il -^ v' ,il') over the scattering angle under the assumption 
hviliv/nii^c^) ^ kTg. In this limit, it turns out to be possible to carry the term describing the effect of Compton recoil out of an 
exponential factor that enters the expression for K(v, fl -^ v' ,il') and implement the subsequent integration analytically. As a result, 
the range of applicability of our approximation for the isotropic kernel is narrower than that of our formula for the angle-dependent 
kernel: hv <^ 50 keV, hv(hv / nif^c^) <^ kT^ <; 25 keV; in this range of parameter values, the accuracy of the approximation is better 
than 98 per cent. The latter constraint means that the average recoil-induced frequency shift must be less than the typical Doppler 
broadening. The opposite case, hv{hv/m^c^) ^ kT^, corresponds to a situation in which the recoil effect is predominant. In this case, 
the singl e-scattering line profile is double-peaked, which is related to the Rayleigh scattering phase function (e.g., P ozdnyakov et 



al. 1979, |1983 ). It is also worth noting that when the temperature of the matter is sufficiently low, scattering of X-rays on neutral 
hydrogen and helium ma y become more important th an Compton scattering on free electrons (see a discussion of the recoil profile 
arising in this problem in Sunyaev & Churazov 1996). 

The paper is organized as follows. In §2 we report our analytical results. For the observational astrophysicist, reading this part of 
the paper should be sufficient for finding all the information necessary for application of the results. The main results are the formulae 
for the K{v,il -> v' ,^') and P{v -^ v') kernels — equations (0),([l^), and equations (|l9l),(0),(||), and (Q, respectively. We 
thoroughly examine the properties of the kernels and determine the parameter ranges of applicability of the different approximations. 
In §2 we also: (1) discuss the properties of the angular function for Compton scattering in a hot plasma that results from our 
K{v,^ -^ v' ,il') kernel; (2) show that the Fokker-Planck e xpansion (to fourth order) of the kinetic equ ation (^ with the kernel 
equation (n9b leads to the generalized Kompaneets equation ( Itohetal. 1998| ; Challinor & Lasenby 1998), equation (W4); (3) derive 
mildly-relativistic formulae, equations ( p9| ) and (51 ), for the kernels, K™°'{i', tl; v' , fJ') and P'™(j/; v'~), that correspond to problems in 
which a decisive role is played by induced Compton scattering; and (4) give the result of the convolution of a spectrum described by 
the step function with P(v -^ v') as an example of application of the kernel to Comptonization problems. 

The subsequent sections of the paper provide the calculation details. The derivation procedure for the K(y, fl -^ v' ,^') kernel is 
described in §3. The angular scattering function is calculated directly [to a better accuracy than K(v^ fj -^ j/', fj')] in §4. The formula 
for the P(y -^ v') kernel is derived by integrating K(v, fJ -^ v' ,il') over the scattering angle in §5. An alternative (fully independent) 
method for deriving the P(v -^ v') kernel in the low-frequency case (hv <C kT^) is presented in §6. 

2. RESULTS 

2.1. The Kii',n^v\n') Kernel 

BR70 have made an attempt to take into account relativistic effects associated with high electron velocities for the kernel 
K(v, n -^ v' ,il'). We have examined the derivation of these authors and found that not all relevant correction terms were included 
by them. In particular, the relativistic corrections to the Maxwellian velocity distribution function were neglected. The BR70 kernel 
is therefore of better accuracy [correct up to terms of order (kTe/meC^y^^hv/m^c^] with respect to photon energy than with respect 
to electron temperature. Therefore, this kernel is strictly only valid for nonrelativistic electrons. In §3, we revise the derivation 
of BR70 and obtain the following approximate formula for KiVjil — > i/'^il'), which is correct up to terms of order (kTg/nieC^)^^-^, 
(kTe/meC^y^^hv/nieC^, and (hiy/meC^)^: 
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(7a) 



where 



g=\vn-v'n'\ = (v^-2vv'n,+v'^y^^ 
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hvv' 
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(7b) 
(7c) 

(7d) 



and Ms = J^fJ'. Note that equation (M gives the probability of a scattering event per unit dimensionless time, t = ajNeCt. 

Let us look at the sequence of terms within the braces in equation (^3). The main term, 1 + Ms' is just the Rayleigh angular scattering 
function. The last term, which is proportional to (hv/meC^y, describes the second-order KJein-Nishina correction to the scattering 
cross section. The remaining five terms become important when the electron velocities are high. These temperature-correction terms 
are either incorrect or absent in the kernel of BR70. The term of order {hv/mf^c^)^ was not given by BR70 either. 



NARROW LINE AFTER SINGLE SCATTERING 



Equation (^ is a good approximation (which will be supported below by a direct comparison with results of numerical calculations) 
to the kernel if both the photon energy and electron temperature are moderately relativistic. Furthermore, this formula without the 
temperature-correction terms (see the resulting eq. [ [lO| ] below) describes the scattering of photons of arbitrary energy (including 
hv ^ m^c^) on nonrelativistic electrons {kTf. ^ m^c^) exactly. 

The term K{v,il -^ v' ,^') must obey the detailed balance principle (including induced effects), i.e., ensure conservation of a 
blackbody spectrum, Bi, = 2hv^ / c\tx^{hv / kTf.)- 1]"', in thermodynamic equilibrium: 
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Equation (M) reduces to 



K{v,^^v' ,^')= [ — I exp 
It is easily verified that our expression (|7|) does satisfy eq. (I 
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2.1.1. K{v, fl —!■;/' , ft') for Nonrelativistic Electrons and Photons of Arbitrary Energy 
If the electrons are nonrelativistic (77 ^ 1), equation (0) simplifies to 
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where g and r] are given by equations ( [7b| ) and ([7d|), respectively. 

Since equation (|l^ fully takes into account the Klein-Nishina scattering cross section, it holds true for photons of arbitrary energy 
when the electrons are nonrelativistic. 

If we are interested in the case in which the photons are of sufficiently low energy, hv <^ 0. Im^c^, the second-order Klein-Nishina 
correction term in equation (|l^ becomes small and can be omitted. The result is 



K^{v,n^ v' ,ii') = --\ -if'i^ —{\+ ^ll)s^^< 
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Equation (njj) includes the first-order Klein-Nishina correction and takes into account Compton recoil; both effects are essentially 
described in the exponential factor (see the derivation in §3). It also accounts for (to first order) the asymmetry of the scattered 
profile due to the Doppler effect (the preexponential factor v' /v in eq. [ 11 1 is important here), and therefore allows one to derive the 
Kompaneets differential equation. 

It is worth noting that the simplified derivation of the K(v, fl -^ v',n') kernel described in an appendix to the BR70 paper makes 
it possible to obtain the exponential factor in equation (jlu), but not the preexponential factor v' /v, and therefore does not lead to the 
Kompaneets equation and does not describe the energy transfer from the electrons to the photons by scattering. 

2.1.2. Angular Scattering Function 

By integrating the K{v,0, -^ v',^') kernel over the photon frequency upon scattering, one can determine how many photons are 
scattered in a unit time through a given angle, i.e., the angular function for scattering (defined by eq. [ JTOJ ] in §4): 



— =27r K{v,n^v\n')dv'. 
d^s J 



(12) 



We have verified, through numerical computation, that the angular function that corresponds to the Km kernel given by equation 
(p^ is exactly the Klein-Nishina formula, describing the cross section for Compton scattering on an electron at rest: 
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Two examples of the angular function that corresponds to the K^ kernel and is described by equation ( |l3| ) are shown in Figure la: 
one for the case hv <C ni^c^ and the other for hv ^ mgC^. We see the well-known Klein-Nishina pattern, namely, that more photons 
are scattered forward (/Xj = 1) than backward (/Xj = -1). This angular function corresponds to the case of nonrelativistic electrons 
(kT, < m,c^). 
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Fig. 1 . — Angular function (in polar coordinates) for Compton scattering on thermal electrons, in various regimes set by the values of parameters: photon energy, 
E = hi/, and electron temperature, T^. (a) Low-temperature case: kT^ <C hu, kT^ <^LmcC^. Solid lines: Angular functions corresponding to the Km angle-dependent 
kernel given by eq. (nOl. These patterns are described by the Klein-Nishina eq. (|l3|). Dotted line: Rayleigh angular function (also in the other panels), (b) Low 
frequency case: hu <^kT^, hu ^m^c- . Solid lines: Results of Monte CarloLsimulations. For the mildly relativistic example, kTc =25 keV, also showjiare the angular 
function calculated from eq. (|l4|), which results from the K kernel (eq. ty) (dashed line), and a more accurate approximation desciibed by eq. ( |l5| ) (dash-doUsd 
line), (c) and (d): Both hu and klc are mildly relativistic. The Monte Carlo results (solid lines) are compared with the results of eq. (|l4) (dashed lines) and eq. (Ila) 
(dash-dotted lines). 



If the electrons are mildly relativistic, the K^- kernel becomes inaccurate, and the temperature corrections included in the K kernel 
(eq. [M]) become important. Analytical integration of this kernel over the scattering angle is possible if an assumption is made that 



hvihv/m^c ) <C kT^ ^ m^c . In this limit, K can be written in the form given by equation ([77|) of §5, which leads to 
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+ 0((hiy/m,c'){kTjm,c'), {kT^/m.cy) 



(14) 



The expression in square brackets in equation ( 14) is an expansion series in powers of hv/nieC^ resulting from equation ( pj| ) (only 
two leading terms are presented). More interestmg is the correction term of order kT^/m^c^. The origin of this term is related to 
Doppler aberration and has nothing to do with quantum effects. The terms that are given in implicit form in equation (14), i.e. 0{...), 
indicate the order of inaccuracy of our approximation for the kernel. Although equation ( p^ was obtained under the assumption 
that hvihv / m^c^) ^ kT^ ^ nigC^, it holds true for arbitrary proportions ofhv and kT^, which we have verified by a direct numerical 
calculation of the integral (p^). 

It is also possible to directly calculate the angular function, not using the Kiy^il -^ v' ,il') kernel. The corresponding derivation 
procedure, which is similar to but simpler than that for K{y^ fl -^ v' ,il'), is described in §4. The result, which is more accurate than 
equation (p4[), is 
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(15) 
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We see that the first-order correction terms, 0{hv / m^c ) and 0{kTi./m^c ), and the second-order term proportional to {hv/m^c ) 
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are the same as in equation (|lj). The last two terms in equation (p3|), 0((hv/meC^)(kTe/meC^)) and 0((kTg/meC^)-^), belong to the 
next-order approximation (with respect to K given by eq. |Q]) to the kernel. 

Integration of equation (|l5|) over all scattering angles leads to the well-known expression that describes the total cross section for 
Compton s cattering on Maxwellian electrons [defined as a = (XNe)~^, where A is the photon mean free path] in the mildly relativistic 
limit (e.g., [Pozdnyakov et al. 1983| [Shestakov et al. 1988|): 



CT = CTt 
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hv 
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(16) 



Note that the pure temperature terms, 0((kTg/meC^)"), in equation (Oh give no contribution to a. This is a well-known fact, which 
means that the total cross section for low-energy photons in a hot plasma is equal to the Thomson cross section. At the same time, 
the angular function (eq. [Oh) is different from the Rayleigh function. 

In Figures lb-Id, several examples of the angular function for Compton scattering on hot electrons are presented; the results of 
Monte Carlo simulations are compared with the results of the calculation by the approximate formulae (n^ and (Ob. 

Figure lb demonstrates scattering of low-frequency photons: hiy <C kTf^, hv ^ nii-c^. The angular function in this case is totally 
different from the Klein-Nishina one — compare with Figure la. Let us first consider the mildly relativistic temperature range 
kT^ <J 25 keV, within which equation (14) (with the terms containing hv/m^c^ vanishing) is a good approximation. Scattering 
is somewhat suppressed (compared to the Rayleigh angular function, which corresponds to the K^ kernel in this case) both in 
the forward and backward directions. There is, however, a noticeable enhancement in the number of photons scattered through 
intermediate angles, between 69° and 138°; these values are found by equating the correction term of order kT^/m^c^ in equation 
(M) to zero. The temperature-relativistic correction to the Rayleigh angular function reaches a maximum of 12(?7/0.05)% at an angle 
of 105°. The relativistic reduction of the angular function is maximal, 10(r//0.05)%, for the two extreme values of the scattering 
angle, and tt. It is evident from Figure lb that the approximation represented by equation ( 15 ) is more accurate that that given by 
equation (|l4). The inclusion of the correction term of order {kT^. / nif-c^)^ is particularly important for very large scattering angles 
(close to tt). Equation (Oh proves to be a good approximation to the actual scattering angular function in the range kT^ <, 35 keV. 

As the temperature becomes significantly relativistic (see the patterns for kT^ = 0.5meC^ and 5m^c^ in Fig. lb), the scattering angular 
function modifies further and becomes totally unlike the Rayleigh angular function. Only the results of Monte Carlo simulations are 
shown for these cases, because the convergence of the expansion series in powers of kTg/nigC^ for the angular function becomes 
poor when kT^^ exceeds ^ 40 keV. Scattering in the forward direction is now heavily suppressed (i.e., the plasma effectively screens 
itself from the radiation incident from outside), while more and more photons are scattered through angles between 7r/2 and tt. In 
particular, at temperatures kTf^ > Q.SmgC^, the number of photons scattered through an angle of tt is higher than in the case kT^^ = 0. 

In the limit kT^ '3> nieC^, the angular function is described by the law da/dfi^ = (1 - /^s)/2 ( Sazonov & Sunyaev 2000[ ). 

Figures Ic and Id show two examples of the angular scattering function when both /iz^ and kT^ are mildly relativistic. The 
approximations described by equations ( [l4| ) and ( [l5| ) work well for photon energies hv <J 50 keV within the temperature ranges 
quoted above. 

The ph enomenon of anisotrop i c (backward) Compton scattering by a hot plasma has been previousl v mentioned i n astrophvsica l 
literature dGhisellini et al. 199l] ; [Titarchuk 1994| ; ^outanen & Svensson 1996| ; [GierHnski et al. 1999| ). In particular, [Haardt (1993| ) 
has performed a semianalytical calculation (using Monte Carlo simulations) of the angular scattering function, as well as other related 
quantities, and obtained results for the case of low-frequency radiation (hv ^ nigC^) that are similar to those depicted in Figure lb. 



2.1.3. The Moments of the Kernel 

Using equation (Ujh of §5 it is possible to calculate the moments of the K{v, il - 
above the angular function. Here we define the moments as 



v',^') kernel, in a similar way as we derived 



{{Av)")f,^ = In jK(v, n -^ v',n')(v'-v)"dv', 



(17) 



so that the integration of {(Av)")^^ over /is gives the moments of the P(v — > v') kernel for the isotropic problem, which were defined 
in equation (Bh. 

The result for the first four moments is 
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The moments of higher degrees turn out to be at least of order Tf . The importance of including relativistic corrections in the kernel 
becomes clear from examining the terms proportional to 77^ and rjhv/m^c^ in the expressions for the first two moments, which 
become comparable to the main terms already at moderate rj^hv/m^c~ ^ 0.02 for large scattering angles. The first and second 
moments describe correspondingly the average frequency increment by scattering and the broadening of the scattered profile. 

It should be noted that the accuracy, to within 0(rf , rjliv/nigC^ , {hv/nii-c^)^), of the derived moments is better than the accuracy, to 
within Oiifl^ ^rj^l^hu /nieC^), of the kernel itself. The reason for this is that K{v, fl — > v' ,il') is multiplied by a small quantity, the 
frequency shift (~ v-q^l^) raised to a certain power, when the moments are calculated. 



2. 1 .4. Comparison of the Analytical Formulae for the Kernel with Numerical Results 

We performed a series of Monte Carlo simulations (a modification of the code described by Pozdnyakov et al. 1983 was used) 
to evaluate the accuracy of our analytical expressions for the K{v,^ -^ v\il') kernel. Figures 2-7 show examples of spectra that 
may form through single scattering (with a given angle) of a monochromatic line on thermal electrons in various scattering regimes 
that are set by the values of the parameters: kT^, hv, /is. The numerical results are compared with the analytical kernels K (mildly 
relativistic) and K^ (non-relativistic). 

As inferred from Figures 2-7 and further from the entire set of results of our numerical calculations, the mildly relativistic equation 
(0) is a very good approximation for electron temperatures kT^ <^ 25 keV and photon energies hv <J 50 keV. In this range of parameter 
values, the accuracy is better than 98%, except in the far wings of the scattered profile. The latter can be roughly defined as the 
regions where |e| ^ 0.5(1 -/is)''^^ (this size should be compared with the characteristic width of the line, which is much smaller: 

[(1 -/is)??]''^^)- For very large scattering angles, /i, <^ -0.8, relativistic corrections are particularly important. In this angular range, 
the accuracy quoted above is achieved in a more narrow temperature range, kTi^ < 15 keV. 

One can safely use the K^^ kernel when kT^ < 5 keV. The photon energy can take arbitrary values in this case because equation (|1( 
is exact for nonrelativistic electrons, as we pointed out before. This is demonstrated by Figures 5-7, which illustrate the scattering 
of a line with hv = 1238 keV (which corresponds to one of the strongest 7 lines produced through radioactive decay of ^^Co during 
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Fig. 2. — Photon number spectra resulting from the single scattering of a monochromatic line of energy, hu = 6.7 keV, on hot electrons with a scattering angle 
of 7r/2, for a number of values ot^jectron temperature. The results of Monte Carlo simulations (solid lines) are compared with various approximations for the 
angle-dependent kernel: Knr (eq. [|lO||, dotted lines) and K (eq. Wn, dashed lines). Note the increasing influence of relativistic corrections (included in the K kernel) 
on the spectrum as the electron temperature becomes higher 
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Fig. 3. — Same as Fig. 2, but for a fixed, mildly relativistic electron temperature, kT,: = 25 keV, and a varying scattering angle. Note that increasing the scattering 
angle for a given temperature has a similar broadening effect on the spectrum as increasing the temperature with the scattering angle fixed (compare with Fig. 2). 
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Fig. 4. — Spectra resulting from the single scattering of a monochromatic line of energy, hu = 6.7 keV, on low-temperature electrons, Tc = 10 K^ for a number 
of scattering angles. The corrections relevant for relativistic electrons are not important in this case. The presented spectra were calculated from eq. fllOl) for the Km 
kernel, which gives the exact result in the case of nonrelativistic electrons and leads to the same results as Monte Carlo simulations. 
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Fig. 5. — Same as Fig. 4, but for high-energy photons, Iw = 1238 keV, and electrons of kT^ = 1 keV. The Xm kernel given by eq. dlCl) is still very accurate at this 
electron temperature and leads to practically identical results with those of Monte Carlo simulations. Note that the total number of scattered photons is drastically 
decreased on going from /i., = 0.5 to ^s = — 1- This is due to the increasing effect of Klein-Nishina corrections on the scattering cross section (compare with the 
nonrelativistic case in Fig. 4). 
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Fig. 7. — Spectra forming through single scattering of monochi'omatic radiation, hu = 50 keV, on mildly relativistic thermal electrons, kT^ = 25 keV, for a sf 
scattering angles. The results of Monte Cai'lo simulations (solid lines) are compared with the different approximations for the angle-dependent kernel: A'm (eq. 
dotted lines) and K (eq. iM, dashed lines). Temperature relativistic corrections, which are included in the K kernel but not in the K^, kernel, are of importanc^n 
this case. 
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supernova explosions) by thermal electrons, for a number of values of the electron temperature and scattering angle. Note that even 
for a temperature of Tg = lO"* K, the Doppler width of the scattered line is ^^ 1 keV for a scattering angle of tt/2. The Ge cryogenic 
detectors of the International Gamma-Ray Astrophysical Observatory (INTEGRAL), scheduled for launch in 2002, will be capable 
of measuring such broadening. 

The line broadening arising from the Doppler effect depends on two parameters: the electron temperature, kT^, and the scattering 
angle, /is- The width of the profile as calculated from equation (jlQ) is strictly proportional to [(1 -/ijjy]'/^. Therefore, in the 
nonrelativistic limit, one can only determine this combination of the two parameters, rather than kT^. and /i, separately, from a 
measurement of line broadening. This point is illustrated in Figures 2 and 3. One can see that identical profiles can be obtained by 
varying either the temperature or the scattering angle. However, if the electrons are mildly relativistic, the correction terms in equation 
(R) become important, with their own dependence on /Xj. This property makes it possible, in principle, to find both kT^^ and fi^ from 
a measurement of a scattered line. However, in many situations, it would be easier to determine /Xj by measuring the recoil-induced 
shift of the line (i.e., the position of the line on the frequency axis, see Figs. 2-7) which is proportional to (1 - iXs)hv/msC^. 

2.2. The P{v —f v') Kernel for the Isotropic Problem 

In the preceeding paragraph we presented the analytical formula for the K{v, fJ -^ v' ,il') kernel, which can be used for describing 
spectra forming as a result of single Compton scattering of a monochromatic line, with a given angle between the direction from 
which photons are supplied (fl) and the observer's direction (fl')- However, in most astrophysical situations we are presented with 
some angular distribution of incident radiation. In such cases, one needs to convolve the K{v^il -^ v' ,il') kernel with this initial 
distribution in order to determine the emergent spectrum. 

If the initial distribution is isotropic, we arrive at the kinetic equation (0), with the kernel derivable by the integration given in 
equation ([ll). This integral can be performed analytically for the K{v,Q, —^ i'' ,fl') kernel given by equation 0n if we assume that 
hv{hv/m^c~) ^ kTg, so that the Doppler broadening prevails over the recoil shift. The details of the integration procedure are given in 
§5 . The essence of this calculation is to write the argument of the exponential in equation (pb as a trinomial. The two terms dependent 
on the photon energy {hv) then turn out to be small with respect to the main term describing the Doppler broadening (because of our 
assumption that hv ^ kT^), and therefore can be carried out of the exponential. The subsequent integration of K over the scattering 
angle becomes straightforward. The result is 
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The kernel given by equation (|l9|) is a series in powers of 7/'/^ (given that /ii^ ^ kT^), written up to the third order 
In the case of scattering of low-frequency radiation in a hot plasma {hv ^ kT^.), equation (p^) simplifies significantly: 
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Equation (EQ), which totally neglects Compton recoil, can be obtained independently from the K(v,il -^ v' ,il') kernel, using a 
simple calculation procedure that employs a transition to the rest frame of the scattering electron (§6). As a matter of fact, we first 
obtained equation ([20|), and this formula allowed us to check (by applying the detailed balance principle; see below) some of the 
terms dependent on hv/nif-c^ in the more general equation (jly). 

We have verified that the expression ( p^ obeys the detailed balance principle, satisfying (in all orders up to rf'/'^) the equation 
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(21) 



which is the analog of the corresponding equation for the K(v, il —^ v' ,n') kernel (eq. [pi]). 



2.2.1. The Kernel Leading to the Kompaneets Equation 

The expression (|l9|) is the sum of four leading terms of the series in powers of r;'/^ for the P{v -^ v') kernel. By retaining a 
smaller number of terms in this series, one can build up cruder approximations to the kernel. Below, we shall consider two such 
approximations. 



The least accurate approximation uses the first term in the series given in equation (19a): 



Po(v^v') = v-'^^7]-'/^po, 



(22) 



with Po given by equation ( |19t. 

The expression ( yJ\ ) is equivalent to the formula obtained by ^^ummer & Mihalas (1967 ). Pq is symmetric in frequency shift: 
PQi~S) = po(S). Therefore, it only describes the Doppler (random-walk) broadening, not accounting for the average increase in the 
photon energy by the Doppler effect. 



A more accurate approximation can be obtained by summing up two leading terms in the series given in equation (19a): 
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(23) 



The Pfi kernel is necessary and sufficient for deriving the Kompaneets equation (g|). For this rea son, we refer to it as the "Kom- 
paneets equation kernel". When hv <C kP^, equation ( |23| ) is equivalent to the formula obtained by Sunyaev (198C ). The Pk kernel 
already accounts for the asymmetry of the scattered line and the corresponding photon heating. The Kompaneets equation kernel 
also takes into account (to first order) Compton recoil. 

2.2.2. The Moments and Normalization of the Kernel 
Important information on the properties of the P(v -^ v') kernel is provided by its moments, which were defined in equation (Bh. 
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The required integration is readily performed for the kernel given by equation ( [1^ after we express the differential dv' appearing 
in equation (||) through dS, 

dv' = 2V2vr]^l\\+2V2ri^l^5 + (i5\ + 2,V25^'q^l^+ ...)d5. (24) 

The quantity 5 was introduced in equation ( |19b| ). The subsequent integration over 5 should be carried out in the limits -oo to oo. 
Collecting terms up to Oijf , rjhiy/mgC^ , (hv /m^c^)^), we obtain 
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(25) 



The mo ments of higher degrees turn out to be at least of order ry' . The values for the moments above reproduce the results of Itoh 
et al. (1998); |ChaIIinor & Lasenby (1998| ) (more general expressions for the first two moments are also available; see Shestakov et 
al. 1988; §4.3 in Nagirner & Poutanen 1994). It is important to mention that the equations (|5]) are valid for arbitrary values of the 
hv/kTe ratio, including the case kT^ = 0, in contrast to the kernel equation (|l9|) itself, which holds in the limit hv{hv/m^c^) <^ kT^ 
only (see the next paragraph for a discussion of the scope of the analytical approximations). 

The moments ( Ga) can also be obtained by integrating over the scattering angle the moments of the angle-dependent K(v,n -^ 
i/',n') kernel (eq. [01). 

The normalization of the kernel equation ( |l9| ) can be calculated similarly to its moments: 
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This expression should be compared with the known expansion series for the total cross section, equation dlN). One can see that the 
terms in square brackets in equation (Eq) describe the inaccuracy of the approximation in equation ( 19 ), whereas the term -2hv/meC^ 



corresponds to the actual first-order Klein-Nishina correction to the cross section (see eq. [|16|]). Note that equation (26h does not 
describe the second-order Klein-Nishina correction to the scattering cross section, although our expression for the angle-dependent 
kernel K(v, il -^ v' ,n') (eq. [R]) does contain the corresponding term. This term was omitted when we integrated Kiv, il —^ v' ,n') 
over the scattering angle (in §5), because in the current limit (hv ^ kT^), inclusion of this term in the formula for P(v -^ v') (eq. 
[u9p) would be inconsistent with the absence of some (unknown) terms of the same order, such as 0(77^) or 0{{hv / m^c^)^ rf^), in this 
formula. 

For the zero-order kernel, Pq (eq. [p2]), the first two moments are 
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The corresponding relations for the Kompaneets equation kernel (E3h are 
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The terms in parentheses in equation (^0|) describe the inaccuracy of the Pk kernel. In particular, the Klein-Nishina reduction of 
the integral cross section is not covered by this approximation. The Pk kernel also neglects altogether the frequency diffusion of 
photons due to Compton recoil (as a result of the recoil-induced frequency shift depending on the scattering angle). The term of 
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order (hiy/mgC^)-^ in the expansion series for ((Ai/)^), which is present in equation ( P5| ) and absent from equation (P9|), is responsible 
for this. The Kompaneets equation (^, which is a Fokker-Planck equation with its coefficients governed by the m oments (P9|), does 
not include the corresp onding dispersion term, the importance of which for the case hi^ '3> kTg was pointed out in ( p^oss et al. 1978 ; 
Illarionov et al. 1979|). 



Having calculated the normalization (g6|) for the P kernel and knowing the exact result for the total scattering cross section (eq. 
[|l6[]), we can try to crudely take into account the terms of order 77^ for P, not calculating them. To this end, we have to renormalize 
the kernel as 

53 f kT,V 19 hv kT, 41 / hv ^ ^"' 



P' = 



10 V WpC^ 



19 hi' kn 41 

5 nigC^ OTpC^ 20 \ rrieC^ 



(31) 



where P isgiven by equation (|l9|). As a result, the normalization of the modified kernel P' is precisely the bracketed expression in 
equation (Q. 

In the case hiy ^ kT^, the kernel given by equation (pO[) can be renormalized similarly: 



Finally, we may introduce the renormalized kernels Pq and P^: 



P^= 


53 / kn Y 

10 \meC^ J 


lels P(^ and P^: 


V 2 m^c^ J 


P'-( 

\ 


5 kT^ Iw 
2 nigC^ jfigC^ 



(32) 



(33) 



(34) 



The assumed normalizations of Pq and P^ are 1 and 1 



2hi'/mgC^, respectively. 
We should stress that the renormalized kernels are still approximations of the same order of uncertainty as the original kernels, but 
these kernels turn out to be more accurate than the original ones, as follows from a comparison with results of numerical calculations, 
which we present in the next paragraph. Note also that the moments are not affected by the renormalization procedure. 

2.2.3. Comparison of the Analytical Approximations for the Kernel with Numerical Results 

The case of \ii> ^ kTe ^ irieC^. In this case, the profile of a Compton-scattered monochromatic line forms through the Doppler 
mechanism alone. The exact kernel can be computed numerically, using equation ( p9[ ) of §6 or by means of Monte Carlo simulations. 
We employ both methods in our analysis. Figures 8 and 9 compare, for two values of the electron temperature {kTi. = 10 keV and 
25 keV), the accurate spectra with the corresponding analytical dependences as calculated in different approximations for the kernel: 
Pt, Po, and Pk (eqs. @,|||3ll). 

As expected, the asymmetry o f the line (domination of the right wing over the left) increases as the temperature grows. One can 
see that the zero approximation ( Hummer & Mihalas 1967 ), Pq, which is symmetric in frequency shift, matches the line profile 
poorly. Therefore, we recommend its usage be restricted to the range kT^ < 500 eV, where Pq is accurate to better than 98%, except 
in the far wings of the line. The latter can be roughly defined as the regions {v' - ly] / (i'' + 1^) ^ 1/4. 

The Kompaneets equation kernel, Pk, which was derived for the considered case (hi^ ^ kTf.) by [Sunyaev (1980 ), works well when 
kTi; <^ 5 keV. At higher temperatures, it becomes important to take into account the relativistic corre ction terms, (1 + \/2Srj^^^)pt, in 
equation (|o|), i.e., to use our most accurate approximation, Pj. At kT^ = 25 keV, the ( Sunyaev 198C ) kernel overestimates the flux at 
the line peak by 10% (see Fig. 9). The flux turns out to be even more overestimated in the wings. The Pt kernel is accurate to better 
than 98% at this temperature. 

At kTg ~ 50 keV, the electrons are already strongly relativistic: ((v/c)^)'/^ ^ 77'/^ ^ 0.5, and the 77'/^ series in equation ( po| ) 
converges poorly, particularly in the wings of the line. At higher temperatures, one must use the exact kernel (which can be calculated 
from eq. [p9|]) and treat the problem numerically. 

Figure 10 is the same as Figure 9, but presents the results for the modified kernels, Pj, Pq, and P^ (eqs. [p2[]-[p4|]). It is seen that 
the renormalization has indeed appreciably improved the accuracy of the approximations. 



Inclusion of quantum effects: hv ^ 0, kTe <C nieC^. The assumption hv{hv / m^^c^) <C kT,^, which we used when deriving equation 
(Il9|), means that the recoil-induced downward shift in the photon frequency must be small compared to the Doppler line broaden- 
ing. Using Monte Carlo simulations, we have established that P and Pk (eqs. [ p^ and [E3|]) remain good approximations up to 
hvihv/nii^c^) ^ kTg, i.e., when the Doppler and recoil effects become comparable. Naturally, the zero-order kernel, Pq, is a very poor 
approximation in this case, because it totally neglects Compton recoil. This is demonstrated by Figures 1 1 and 12, which compare 
accurate (Monte Carlo) single-scattering spectra with the corresponding profiles computed from the renormalized kernels P', Pq, and 
P'k (eqs. J^, [||, and Q). 

The results of our simulations imply that the temperature limits quoted above for the case /ii/ <C kT^^ remain valid upon inclusion 
of quantum effects. We finally conclude that the P kernel can be used in the following range of parameter values: hvihv / m^c^) <^ 

kTf^ <; 25 keV, hv <^ 50 keV. The corresponding limits for the Pk kernel are hvQiv / m^c^) <^ kT^ <^ 5 keV, hv <^ 50 keV. 
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Fig. 8. — (a) Spectnim that forms through single scattering of isotropic monochromatic low-frequency (hv <C /iTi) radiation on weakly relativistic theoual 
electrons, kT,^ = 10 keV. The soHd Hne showithe result of an accurate numerical calculation from eq. (p9|).^so shown are the approximations given by: eq. (E2), 
the zero-order kernel Pq (dotted line); eq. (ES), the Kompaneets equation kernel Pk (dashed line); and eq. fctj), the mildly relativistic kernel Pj- (dash-dotted hne). 
(b) Ratio of the approximate spectra shown m (a) to the accurate spectrum. 
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Fig. 9. — Same as Fig. 8, but for a higher temperature, kT^ — 25 keV. 
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Fig. 10. — Same as Fig. 9, but the renormalized kernels P^ (eq. [B3p), P^ (eq. |p4[), and Pj (eq. ||32|]), are used. One can see that the agreement between the 
approximations and the exact kernel is better than in Fig. 9. 
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Fig. 1 1 . — (a) Spectrum that forms through single scattering of isotropic monochromatic radiation of high energy, hv = 50 keV, on weakly relativistic thermal 
electrons kT, = 10 keV. The result of a Monte Carlo simulation (solid line) is compared with the different approximations for the kernel: P^ (eq. [R3||, dotted line), 
P'l, (eq. |R4p, dashed line), and P' (eq. fpm, dasli-dotted line), (b) Ratio of the approximate spectra shown in (a) to the accurate (Monte Carlo) spectrum. 
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Fig. 12. — Same as Fig. 1 1 {hv = 50 keV), but for a higher electron temperature, kT^ = 25 keV. 



At hv{hv/m^c^) ^ kT^, recoil becomes more important tha n the Doppler effect. The single-scattered profile becomes double- 
peaked (see, e.g., the results of Monte Carlo simulations of Pozdnyakov et al. 1983). This case requires a special analytical 
treatment that is beyond the scope of this paper Here we may only point out the principle mathematical diffuculty that makes it 
impossible to use our current method for calculating the P{v -^ v') kernel in this limit; the exponential entering the expression for K 
(eq. [R]) cannot be expanded in powers of 77'/^, as was done for the case hvihv / m^c^) <^ kT^. 

Numerous examples of spectra forming by single Compton scattering, which may well be encountered in various astrophysical 
situations, and which are well approximated by our analytical kernel equation (Jl9|), are presented in Figures 13-15. 

2.2.4. Properties of the Single-Scattering Profile 

The line profile forming by single scattering of isotropic radiation on thermal electrons (see Figs. 8-15), which is approximated 
by equation (p^, is unique in its properties. It therefore appears interesting to examine its characteristic features in some detail 
(following Pozdnyakov et al. 1983). 

First, let us compare the line profile as calculated from equation ( jl^ ) with the usual Gaussian profile that, e.g., may result from 
Doppler broadening of an emission line in the presence of thermal or turbulent motions of ions. To simplify this comparison, 
let us assume that hu <C kT^^. For a given plasma temperature, it is natural to adopt /S.vo = i/(2?7)''^ for the Doppler shift[], i.e. 
A^'^ exp[-(j/'-i/)^/Af||]. The mean (rms) frequency shift, ^((Aj/)^), is lyrj^^^ for the Doppler profile. The corresponding value for 
the Compton-scattered line is i/[277(l +23.5r7)]'/^ (as results from the value of the second moment given by eq. [p5|]). The width of 
the single-scattering profile at half-maximum (FWHM) is approximately equal to Ivllnlrj]^^^ = 1 .661^77'''^. The corresponding value 
for the Doppler profile is 2v[2\\\2ri\^l^, i.e., V2 times more, which is opposite to the situation with the rms shift. Thus, in the case of 
the Une forming by Compton scattering, relatively few photons appear in the upper part of the profile (above half-maximum), and an 
accordingly large fraction of the scattered radiation emerges in the wings of the line. It is also worth noting that the Doppler profile 
is symmetric, while the profile due to Compton scattering is not. 

Now let us consider the peak of the single-scattering profile, a detail which makes it so peculiar. In the vicinity of the maximum 
(|i/'- z^l ^ vq^'^), the spectrum is well approximated by the following expression, which results from equation (^ 



P(v -^ v')+- = V 



11 

20 V TT 

v' lv-\ 



■1/2, 



1 + 



1091 23 



616 154 V kT^ 



V 



1/2 



22 V 2 2 \ kT'' 




(35) 



'Note that the width u.sed, Aud, is (M/nic)'^^ = 43(M/mp)''^ times the thermal width of hnes of an ion of mass M. 
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Fig. 13. — Spectra resulting from the single scattering of isotropic monochromatic radiation of energy hv = 6.7 keV on low-temperature {hv > 4kTs) thermal 
electrons, for different values of kT^ . In this case, the Compton-recoil ,sliift is larger than the Doppler shift. The results of Monte Carlo simulations (solid lines) are 
compared with the results of the calculation by the approximate eq. (Hll) for the mildly relativistic kernel P' (dashed lines). For the case kT^ = (cold electrons), 
only the Monte Carlo result (double-peaked profile) is shown, since our approximation for the kernel is not valid in this case. 
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Fig. 14. — Same as Fig. 13 (hi/ = 6.7 keV), but for high-temperature electrons, hu < 4kTe. In this case, the Doppler shift is larger than the recoil shift. 
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Fig. 15. — Spectra resulting from the single scattering of isotropic monochromatic radiation on weakly relativistic electrons, kT^ =jji) keV, for different photon 
energies. The results of Monte Carlo simulations (.solid lines) are compared with the results of the calculation by the approximate eq. (pll) for the mildly relativistic 
kernel P' (dashed lines). One can observe how the effect of Compton recoil on the spectrum increases as the photon energy becomes higKer. The case of the 122 keV 
nuclear line produced by ""^Co is beyond the scope of our analytical approximation for the isotropic kernel. 

where the indices plu and minus signs correspond to the right and left wings, respectively. 

We see that the spectrum has a cusp ?Av' = v (a break in the derivative occurs there). Near the cusp, on its both sides, the spectrum 
can be approximated as a power law, the slopes in the right and left wings [the coefficient at (v' jv— 1) in eq. (p5|)] being significantly 
different: 



a+ 



dXviP 



d\xv{v'lv) 
d\a.P 

d\n{v' jv) 



=1^+0 



v'=v-% 



22 V 2 ' 2 2kV 



22 V 2 ' 2 2kT' 



-0L+-- 



hv 



(36) 



(in Pozdnyakov et al. 1983, a_ - a+ = 3 when hv = 0, because they considered the energy spectrum that is the product of ly' /ly and 
the photon spectrum considered here). 

It is interesting that when hi/ = kT^, the line profile in the vicinity of the cusp is symmetric (in logarithmic coordinates) about v' = 1/ 
ia+ = a-). Let us give a few further examples. lfhv = Q and ?/ = 0.001, 0.01, and 0.1, the slopes run: a+- = 27.0 =F 0.5, 8.5^0.5, and 
2.7 T 0.5. For hv/m^c^ = 0.1 and 77 = 0.02 (the case of Figure 1 1): a+ = 8.0, a. = 4.1, and a. - a+ = -3.9 (for /ji^ = and 7/ = 0.02, 
we would have a+ = 5.5, a- = 6.6, and a--a+ = 1.0). 

The asymmetry of the single-scattering profile can be demonstrated by comparing the fractions of the scattered radiation contained 
in the right {v' > u) and left {v' < v) wings of the line. Using the renormalized kernel given by equations (BTI) and ( 19 ), we find in 
terms of number of photons (N = J P' dv') and total energy (W = / hv'P' dv'). 
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(38) 



The renormalization is important here because it enables us to obtain the terms 0(1]^), 0{rjhv/meC^), and 0((hi'/meC^)^) in equations 
(BtI) and (Esb. This procedure is strictly correct for the following reason. The terms of even orders in 77'/^ in the power series for 
the P(t/ ^ I'') kernel (see eq. [p^), i.e. po, pt, p,- and analogous (unknown) terms of higher orders, are symmetric in frequency 
variation. Therefore, if we know (and we indeed do) the contribution of such a term, say 0{rf), to the normalization of the accurate 
kernel, we immediately know that the contribution of this term to both A'+ and A'- is equal to half this value. W+- are then determined 
as ^ h{v + v' - v)P' dv' = N+-hv + j h{v' - v)?' dv' . The last integral is accurate to within rf-, rjhv/m^c^, and {hv/m^c^)^ because of 
the presence of the small factor v' - v. 

An interesting conclusion can be drawn from both equation (B8h and the above discussion: the left wing contributes to the total 
accumulation of energy by the photons exactly as much as the right wing. 



2.2.5. Single Scattering of a Step-Function Spectrum 

Although the present paper is mainly devoted to the case of single Compton scattering of narrow spectral lines, our basic formulae 
— equation (^ for the direction-dependent problem and equation ( p^ for the isotropic problem — describe the kernels of the 
corresponding integral kinetic equations (eqs. '^ and |^]) appearing in the general problem of Comptonization in thermal plasmas. 
Let us present a simple example of using the P(v -^ v') kernel. Consider the scattering of an isotropic photon distribution described 
by the step function, i.e., dN^/dv = 1 if t^ < i/q and dN^/dv = if t^ > vq, in an optically thin (so that multiple scatterings are not 
important), hot plasma. The spectrum of the scattered photons can be found by convolving the initial frequency distribution with the 

dN{v) fdN^iv') , 

-P^v —f u)di' , (39) 



du 



dv 



where r ^ 1 is the Thomson optical depth of the scattering medium. 

Figure 16 presents examples of spectra of the scattered component as calculated from equation ( p^ using different analytical 
approximations for the kernel: Pq (eq. [p3|]), f^ (eq. [p4|]), and P' (eq. [pi|]). The radiation is assumed to be low-frequency 
{hi'() <C kTs). In one case, kT^ = 1 keV, a spectrum forming in the case of non-negligible photon energy {hv^ = 7.1 keV) is also shown. 
The integral in equation (B9h was performed numerically. Note that the actual observable spectrum will be the sum of the unscattered 
and scattered components, and will depend on r, i.e., {dN / dv\otii\ = {\-T)dNi)/dv + TdN /dv. 

The spectra shown in Figure 1 6 allow analytical description. We present here only the result of the convolution of the step function 
with the Kompaneets equation kernel, Pk- 

dv \ (1-/o + /i)t, v<vq, 



(40a) 
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(40b) 



The main t erm in equation (|40a|) — /qt (if t^ > t^p) or (1 - /'o)r (if v < t/p) — results from the zero-order kernel, Pq ( Hummer & 
Mihalas 1967). 

The spectrum of the scattered component at frequencies t^ > I'p has a quasi-power-law shape, with a slope that is approximately 
equal to the slope of the right wing of the kernel itself (see eq. Pq]): 
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(41) 



As is the case with the kernel itself, the first-order temperature correction causes the spectrum of the scattered component to be flatter 
than results from the zero-order approximation. Compton recoil has an opposite effect on the slope and may cause a significant 
steepening of the spectrum when /zi'o ^ kT^. It is clear from Figure 16 and equation (^Tj) that in the low-frequency case, it is possible 
to determine the temperature of the scattering plasma by measuring the slope of the spectrum of the scattered component. 

Let us also give an expression for the total number of photons that have been redistributed from the frequency region v <vq into 
the region v > V(j: 
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(42) 
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Fig . 16. — Single scattering of a specttum described by the step function — iINq jdE = 1 if £ < £o and dA'o jdE = 0i{ E > Eq — on mildly relativistic themial 
electrons, considered for various electron temperatures, assuming that the photon energy is negligible (Eq <IC kTi:). The spectra pesulting from the convolution (eq. 
|B9||) of the step function with the following approximations for the isotropic kernel are shown: the zero-order kernel P^ (eq. |p3p, dotted lines), the Kompaneets 
equation kernel P^ (eq. [R4], dashed lines), and the mildly relativistic kernel P' (eq. |pl|], dash-dotted lines). For the case of kTc = 1 keV, also shown is the spectrum 
(corresponding to the P'kernel) for a nonnegligible photon energy, Eq = 1.1 keV (solid line). Note that the actual observable spectrum will be the sum of the 
nonscattered and scattered components, i.e., (dN/dE)fau\ = (l—T)dNQ/dE + TdN/dE, where r <C 1 is the Thomson optical depth of the scattering medium. 



We can also imagine a situation in which a spectrum described by the left-side step function: dNo/di' = lif u>i^o and dN{)/dv = 
\fv< Vd, is being scattered. In this case, we arrive at a formula that is similar to equation (40): 



^^ f (./o-./i)t, y<ui.) 
dv 1 (1-./o-/i)t, v>vq. 



(43) 



We see that here the asymmetry of the kernel has the opposite effect on the scattered spectrum as c omp ared with the previous 
situation, in which the right-side step function was considered (compare the signs of the term f\ in eqs. [40a] and [P3|]), namely, the 
temperature correction causes the spectrum to steepen, while the recoil correction makes the slope flatter. 



2.3. Fokker-Planck Expansion of the Integral Kinetic Equation 

We can carry out a Fokker-Planck expansion (eq. [0]) of the kinetic equation (S) with the P(i' -^ u') kernel derived in this paper 
(eq. [p^). The coefficients in equation (Eh depend on the moments of the kernel, which are given by equation (ESb. As a result, we 
obtainthe generalized (for the mildly relativistic case) Kompaneets equation: 
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(44) 



This equation was earlier obtained, in a different way, by Itoh et al. (1998) and Challinor & Lasenby (1998 ). Substituting the 
Planckian distribution that corresponds to the temperature of the electrons, n = (e^"^/*' ° - 1)"', into equation (|44|) yields dn(v)/dT = 0. 
This test confirms once more that equation ( p^ ) for the kernel consistently takes into account necessary corrections. One can also 
make sure that equation (H) without the terms responsible for induced scattering does not modify a Wien spectrum, n = 6"''''/*^= . The 
fact that equation (Ej) conserves the total number of photons follows directly from its divergent form. Thus, as expected, the basic 
properties of the Kompaneets equation are retained in equation (p4|). 
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Equation ( p4| ) can be simplified in the case hiy ^ kTg. Omitting the terms related to quantum effects and induced scattering, we 
find 
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(45) 



This equation allows one to derive the first-order relativistic corrections to the effect of distortion of the CMB spectrum in clusters 
of galaxies with hot gas ([toh et al. 1998; phallinor & Lasenby 1998;see another way of finding these analytical corrections in 



Sazonov & Sunyaev 1998)^ 

In another limit, hvihv / mgC^) ^ kT^, when the Doppler effect is not important, one obtains (ignoring induced-scattering terms) 
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The second parenthesized te rm i n this equation, which describes frequency diffusion of photons, was added to the Kompaneets 
equation by Rossetal. (1978)and Qlarionovet al. (1979). This correction becomes especially important when studying the scattering 
of hard radiation on cold electrons in an optically thick medium. Such a situation takes place, for example, during a supernova 
explosion; an analytical solution to the corresp onding diffusion problem wa s derived and employed in a calculation of the evolution 
of the X-ray spectrum of Supernova 1987A by Grebenev & Sunyaev (1987| ). 

In the intermediate case of hv ^ kT^., hv{hv / nigC'^) ^ kT^., the dispersion term (kTg/mgC^)v~^d/dv{v'^dn/dv) of the Kompaneets 
equation, which describes the diffusion of the photons due to the Doppler effect, must be added to equation (p5). 

2.4. Kinetic Equation for Problems with a Decisive Role of Induced Compton Scattering 
If the conditions n^ 1 and n ^ kT^/hv are both satisfied, then equation (Q) will simplify to (only induced-scattering terms have 



been retained) 
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It is interesting that only correction terms that are proportional to kTg/mgC^ appear in this equation. There is no term proportional to 
hv/mgC^, although such a term is present in the diffusion equation describing the spontaneous scattering process, equation (^6|). This 
is a result of the joint operation of Compton recoil and Klein-Nishina corrections, which both contribute to the first two moments of 
the kernel (eq. [g|]). 

In the past, many phenomena caused by induced Compton scattering were investigated in the nonrelativistic approximation, using 
only the main te rm in equation (p7|). One such phenomenon is distortions in the low-frequency radiation spectra of radio sources 
(Sunyaev 1971), which become large if kTt, = O.S/^A^ ^ meC^/r(l +t), where /^ is the intensity of quasi-isotropic radiation at a 



wavelength A and t is the Thomson optical depth of the scattering cloud. Particularly interesting is the case of bright extragalactic 
radio sources, for which kT/, ^ nigC^ /t eve n though r <C 1, because T i , ^ 10" -^ 10'^ K. Other phenom ena include plasma heating 
( Peyraud 1968|;|Zeldovich & Levich 1970|; t^evich & Sunyaev 1971 ; Vinogradov & Pustovalov 1972|; B landford & Scharlemann 
1976; plarionov & Kompaneets 1977| ) and induced light-pressure force ( |Levich et al. 1972| ) in the vicinity of astrophysical objects 
emitting low-frequency radiation with high Tt,. Obviously, relativistic corrections, described by equation (|4^), could play an important 
role in such problems. In particular, these terms (although small) should play the ro le of viscosity for such phen omena as the 
formation of shock waves in the photon spectrum during Bose-condensation of photons ( Zeldovich & Sunyaev 1972 ). 

Induced Compton scattering may also lead to essentially aniso tropic effects, such as narrowing or spreading ( depending on the 
spectrumof the radiation) of a radiation beam traversing a plasma (Goldinet al. 1975; Zeldovich & Sunyaev 1976| ). It may also play 
a major role in the interaction of beams of maser radiation having narrow spectra (Av < Avo, where Avo = t^[2(l-/is)fc7'e/'«eC^]''^ 
is the Doppler shift) (Zeldovich et al. 1972). Let us write (for an infinite homogeneous medium) the integral kinetic equation that 
arises in such problems. 



dn{v,il,T) 
d^ 



n(v,n,T) n(v',il',T) 



K(v',n' -^ v,n)-K(v,n'^v',n') 



dv'dil' 



= n(v,n,T) / n(v' ,n' ,T)K''"^(v,n;v' ,n')dv'dn' , 



(48) 



where n(v) is the occupation number in photon phase space, and we have introduced the new kernel K"'^{v,n;v' ,il'). In the limit 
hv ^ ififC^, hv ^ kTe, which always takes place in the case of compact sources of low-frequency radiation, a particularly simple 
expression for this kernel can be given. Indeed, we can expand the exponential in the expression ( [7a| ) according to equation (p3) of 
§5, and then, taking the second term in the resulting series out of the exponential, get 
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Equation (S^, like the formula in equation (^ for the K{v^ fl -^ v' ,^') kernel, can be used when the plasma is mildly relativistic. 



kTf. <J O.lmeirT In the case of nonrelativistic electrons, kT^. ^ mgC , the correction terms in equation ( [4-9[ ) can be neglected, and 
can be replaced by v everywhere except in the difference v' - v, which results in the following formula 

{v'-vf 



K^\u,n-u',n')=-{T:r,) 



lr^.rV2,.-iM^ 



rrigC^ 



(l-/xsr'/2(l + //^)exp 



4:^Hl-tJ-s)v 



(50) 



which was earlier obtained byl Zeldovich et al. (1972[ ). 

Integrating the kernel (eq. [[4-9|]) over all scattering angles gives the corresponding kernel for the isotropic problem. This kernel is, 
however, more easily deduced from equation (U^: 



p'^d^j^.j,')^ 



P{v' -fv)-P{v~^v') = 



^ 2 2/ _, A (PO + P')i 



(51) 



where po and p, are given by equation (19b). 

The kernel given in equation (51) is applicable in the limit hv <C kT^ <J O.lmgC^ and makes possible to obtain the differential 
equation (EtI). 

3. DERIVATION OF THE K(u, fJ -^ ly'.il') KERNEL 

Consider a photon of frequency i' that propagates in the direction Q, . We calculate the probability (per unit dimensionless time, r) 
that the photon will be scattered by a Maxwellian distribution of electrons into a solid-angle interval d^', with the emergent photon 
frequency falling in an interval diy' . 

Our primary goal is to derive a formula that would be a good approximation for situations in which both the electrons and photons 
are mildly relativistic, i.e., when 77 = kT^ / m^c^ , hi' / irieC^ ~ O.L Therefore, we make an initial assumption that rj^hv/mg^c^ <C m^c^. 
The final formula will contain correction terms of order up to 77^/^, rj^^^lu'/meC^, and {hv/m^c^f-. The term of order {hv/m^c^)^ 
originates directly from the Klein-Nishina formula for the scattering cross section. We retain this term in our final expression 
on purpose (although it causes the formula to be of a slighly better accuracy in terms of photon energy that in terms of electron 
temperature), because this expression without the temperature-correction terms then becomes the exact formula for the case of 
scattering of photons of arbitrary (including hv ^ mgC^) energy on nonrelativistic electrons (?7 <C 1). 

Let us introduce the following system of reference: Ox points along fJ, Oy is in the ( D, il') plane, Oz is norm al to this scattering 
plane. There are two basic equations. The first is the energy-conservation relation (e.g., Pozdnyakov et al. 1983): 
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(52) 



where c/3 is the electron velocity and cos a = ilil'. The second equation describes the differential cross section for Compton scattering 
([Fauch & Rohriich 19761; iBeresteskii et al. 19821): 
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Introducing the components of the electron velocity {/3x, /3y, /Sj), we find 

ilj3 = /3.V, n'/3 = /3a- cos a + (3y sin a. 
Equation ( p4 ) imposes a link between the different /3 components for given v' /v and a: 
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In order to calculate K(v, fl —■ 

Kiv,n 



l-(/?2 + /32 + ^2)- 

v' ,il'), we ought to carry out the following integration over electron velocities: 
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(53) 



(54) 



(55) 



(56) 



(57) 



(58) 



Here, the factor (1 - 0/3) allows for the relative velocity of the photon and electron along the direction of the former's motion (L andau 
& Lifshitz 1975), f0{/3x,f3y,l3z) is the electron velocity distribution function. In equation (p8[), one of the velocity components (say. 



24 



SAZONOV & SUNYAEV 



f3y) must be expressed through the other two. The interval, di/', of the photon frequency after scattering can then be related with the 
corresponding interval of Py-, hence the appearence of the factor |9/3j,/9i/'| in equation (pq). 
For fp we substitute the relativistic Maxwellian distribution function, 



//3 = (27TT]) 



-3/2 



1 + — 77 
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3^ 
87 



exp 



277 



(59) 



Here we have retained only correction terms of order r/. 

We can proceed with the integration in equation (Em on expanding the integrand in powers of /3. The correct account of temperature 
terms of order /3^ and P^ necessitates inclusion of the corresponding terms in equation (p9b for the velocity distribution (contrary to 
the statement made in BR70). We also note that the terms of order f3^, (3hv/msC^, and {hv / mgC^Y , which we are keep throughout, 
were neglected altogether in the derivation of BR70. 

As follows from equation (pq), the derivative dj3y/dv' to the first order is 






-d-A). 
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The last bracketed term in the expression 
1 /7. Using equation (pi) we find 



which finally yields 



gives rise to additional terms of order jShv/m^c^ due to the presence of the factor 

5(1/7) 
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(62) 



Note that the main Klein-Nishina correction to the scattering cross section, which is of the order of hv/m^c^, is contained in the 
factor {v' jvf- in equation ( p3h , rather than in the X function (eq. |54|), which describes Doppler aberration and the Klein-Nishina 
correction of order {hv/mJrY'- The reciprocal factor, v/v''^, enters equation (62 ), so upon multiplication of da/dn' and \df3y/di''\ 
in the integrand of equation (pq), the presence of the Oihv/mf.c^) correction in Kiv, fl -^ v' .^') is not explicit. 

Proceeding further with expansions, we obtain 



(Tj dll 



dPy 



dv' 



— : — {\+^ll+ 2s(-/is + ^ll) + 2/3^(-/is + ^ll) 

loTT 1/sma 

+i^(l-Vs + 3Ais) + 2iA(l-3M.s + 2/is) + /3x(l-Vs + 3A*s) + /3Vl+2Ais-3Ms) 



+2/32i(-l +4Ai,-3/i2) + 2/32/3,(-l + Vs-3a*s)+ ■ 



hv 



hv 



r(l-/^s) 



hv n-jx^KUj^ v' ( hv 



(1 + M^) 
1+Ms 



(l-A^s)' 



(63) 



where /is = cos a. 
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Here P^ must be substituted for (3^ in equation (p3h and in the terms preceeding the exponential in equation (p9b. 

Note that in the situation of scattering of nonrelativistic electrons (r; <C 1), equation (tel) reduces to a simple formula: 
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(66) 



which is valid for arbitrary values of photon energy, including hv 3> nieC^. Our final formula will consequently be exact in the limit 
77 ^ 1, as we claimed at the beginning of this section. 

Having made this remark, let us return to considering our main situation of interest, i.e,. with mildly relativistic electrons and 
photons. The accuracy of equation (p5|) turns out to be insufficent for describing the exponential factor of the distribution function 
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(eq. [p9[|), in which (3^ is divided by ry. In this factor, 0^ need be given accurately to the fifth order, which requires the inclusion of 
the factor I/7 in equation ( pq ) (similarly to the situation with dPy/dv' above). We consequently get 

, ^ , , 4^+A^\ / A/32 
^''^P 1 "— ) = exp I — — I ~ ( 1 - -^rr ) sxp I "— I ' (67) 

where 




If] I \ If] 




1— cosa hv 
sin a nieC^ 



A/3^ = -^^^^^/3,/3^ (68) 



with /3y and f]^ given by equations ( pq ) and (p5h, respectively. 

Now, having completed all the necessary preparations, we can implement the integration in equation (B^). Integrals of the following 
type then appear: 



where a and b are constants set by equation (p5h. Such integrals are readily done (see, e.g., BR70). It is natural to present the final 
result in the form of a series in terms of the quantity 

A = i/'-i/+zy'^^(l-cosa). (69) 

OTeC 

Indeed, consider the situation with kTf^ = 0, when all scattered photons undergo the same decrement in frequency, owing to Compton 
recoil: i^/i^' = 1 + (1 - cos a) hiy/meC^, as follows from equation (p^). This shift corresponds exactly to A = 0. If we now allow a 
nonvanishing electron temperature {kTf. ^ 0), the scattered profile will be Doppler-broadened near the recoil-shifted peak of the line. 
The term A will then measure the frequency variation relative to this peak and therefore will always be of the order of vr\^l'^, which 
ensures that the final expression will be convergent regardless of the proportion of hv and kT^ (for comparison, we can consider a 
similar quantity v' - v, which is ~ 1/77'/^ if the Doppler effect is dominant but ^ -hv^ Inyc^ if recoil is more important). Moreover, 
note that A enters as an entity in all the major expressions we obtained above (see eqs. [^], [|62|]-[p5|]). After lengthy calculations, 
we finally arrive at equation (0) given in §2. 

4. CALCULATION OF THE ANGULAR SCATTERING FUNCTION 

The angular function, da/dfi^, for Compton scattering on Maxwellian electrons can be directly calculated using the same for- 
malism that we employed to derive the K{v,il — > v\^') kernel in the preceeding section. In fact, the derivation of da/dfi^ is less 
tedious than that of K(iy, il -^ i^' Al'). For this reason, we here calculate the angular function to a better accuracy than that of our 
final formula for the kernel (eq. [^]). Namely, we are willing to obtain an expression that will contain correction terms of the order 
of (kTg / nigC^)^ and {hv / m^c^){kTs / m^c^) to the Rayleigh angular function. The term of order {hv/nisC^Y will also be found, but this 
term already follows from equation (jlQ) for the K^^ kernel, which, as we recall, is accurate for arbitrary photon energies in the case 
of nonrelativistic electrons {kT^ ^ rn^). 

In order to find the angular function, the integration over the electron velocity space needs to be done: 

^ = ^ y ^(1 - n/3)//3(A, Py,[iddMMP,. (70) 

Here the differential scattering crosscsection, da/dVl', for given ((3^,Py,(3v) and hu/nieC^ is given by equations ( p2[ ) and (p3|), and 
the relativistic Maxwellian distribution function is represented by the series given in equation (p^). The principal difference (which 
simplifies the calculation) of equation (pQ) from the similar equation (Bq) is that /3,, is now a free parameter, like the other components 
of the electron velocity. 

Next, the quantity {\-il(3)d(j/dil' needs to be expanded in powers of (3x, /S,., /3- (to fourth order) and hv/m^c^ (to second order), 
as we similarly did (to a worse accuracy) in the previous section (see eq. |p3]). The resultant expression is rather cumbersome, so we 
do not give it here. For the fp distribution function, the approximation of equation (|5^) is sufficient, because the next order terms, 
0{rf ,ri0^ , (3'^ , (3^ / rj , (3^ / rf), in the series fp cancel upon integration. 

The integration in equation ( [70| ) is connected to the calculation of standard integrals 

P\f3\,i3':^ exp I - ^" ^2 dMPydPz, (71) 

where k, I and m are even numbers [odd terms with respect to one of the (/3f ,/?,., /S^) components vanish upon integration]. 



The final result is equation ( 15 ) 
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5. DERIVATION OF THE P(!^-^!^') KERNEL 



m,. 



Here we demonstrate how to perform the integral in equation (|l]) with Kiv^ f! -^ v' ^^') given by equation (0) when hv{hv/meC^) <C 



Write the exponential factor entering equation ( [7q ) in the form of a polynomial: 

e^ (j/'-i/)2 hvv' {v'-v){l-n,) fhi^i^'Yil-^j,,) 
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(72) 



where 

g2 = 2z.i.'(l-^,) + (z.'-i/)2. (73) 

One can see from equation ( [73[ ) that g-^ ^ at v' = v when /is ^ 1 . The factor g^ enters the denominator of the first member of the 
polynomial in equation (^2b, which describes the Doppler broadening. Knowing this property, which means that K(iy, il -^ v' ,^') is 
a 5-function for /is = 0, one can immediately make a prediction that the P(i/ -^ v') kernel resulting from the integration in equation 
(|l]) over /is will have a cusp (a point where a break in the derivative occurs) at v' = v. Such a cusp is indeed present in our final 
expression for P(v -^ v'), as we shall see below. 

There is no singularity at /ts = 1 in the second and third members of the polynomial ([74), which describe the frequency variation 
due to Compton recoil. This suggests that the cusp mentioned above will remain at the same position, v' = v, regardless of the 
initial photon energy ihv ), which indeed proves to be the case (see the results of Monte Carlo simulations in Fig. 1 of the review by 



Pozdnyakov et al. 1983). 



By assumption {hv ^ kTf^), the main contribution to the photon frequency increment by scattering comes from the Doppler effect. 
This implies that v' -v \^ typically ^ vq^l^. The second and third terms in equation ( [72| ) are thus infinitesimal of the order of 77'/^ 
and ?7, respectively. One can therefore take these terms out of the exponential, which will make possible analytical integration in 
equation (|]). 

It is convenient to present our final result for K{v, il—^v', Jl') in terms of the quantity 



^ = (2?7)-i/2 



v+v 



(74) 



This new variable, describing the relative frequency shift, is similar to e, in powers of which the original equation (R) is written, 
when hv -^ (see the main term of eq. [ f7^ ] below). Moreover, the combination {v' -i')/{i'+i'') arises in a natural way if one first 
calculates the analog of the P{v -^ v') kernel for a monoenergetic isotropic electron distribution and then convolves this quantity with 
the Maxwellian distribution function (as we do in §6). Indeed, for a given electron speed, j3c, the photon frequency after scattering 
can take values in the range \v' -v\/(v+v') < (3. From this fact it becomes immediately clear that P(j/ —> v') must (and indeed does) 
possess a certain symmetry if expressed in terms of 5 given by equation (^3). 
After a few intermediate steps (see eq. [^), which are 
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we finally obtain 
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Note that we have omitted in equation ( |77| ) the second-order, 0{{hv/m^c^)^), Klein-Nishina correction term, which was present in 
the original expression in equation (^. This is because in the limit we are currently working with {hv ^ kTe), inclusion of this term 
would be inconsistent with the absence of (unknown) terms of the same order, such as 0(ri^) or Oiihv/m^c^f'rf^ ), in the series given 
in equation ( |77a| ). 

The term K{v,il 



v\il') in the form (77) is easily integrated over the scattering angle using the change of variables (/is 
2(5^/(1 -/j..s) = t^, which results in the final equation ( [l^ for the P{v -^ v') kernel for the isotropic problem. 

6. DIRECT CALCULATION OF THE P{u -^ u') KERNEL IN THE CASE OF Hu < KT^ 



t) 



Given is an isotropic field of electromagnetic radiation of frequency v. Its spectrum (number of photons per unit solid angle, unit 
frequency interval, and unit detector area) is 

drtdv A-ir 

Consider scattering of the radiation on an electron moving at a speed of v = (3c. The energy of the photons is assumed to be 

low enough (hi' <C WeV^) that Compton recoil can be ignored. We consider electrons that are not too relativistic: (hi>/meC^)j ^ 1, 

where 7 = (1-/3^)"''^; thus Klein-Nishina relativistic corrections are not significant. For the moment, we ignore induced scattering. 

Consider the situation in the electron rest frame. In this frame, the spectral intensity of the incident radiation is direction dependent: 
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(79) 



where /i is the cosine of the angle between the velocity of the electron and the direction of propagation of the photon. Quantities that 
are measured in the electron rest frame (with subscript "0") are related to the corresponding quantities measured in the laboratory 
frame via the Lorentz transformations: 
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Under the given assumptions, the scattering can be treated in the Thomson limit (the photon frequency does not change) in the 
electron rest frame. Therefore, the number of photons scattered into an interval dH-o of solid angle in a unit time is 
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The reverse transition to the laboratory frame can be performed using equation (^ 
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(83) 



where u = v/v. 

The derivation given above is analogous to the one in ( [Sazonov & Sunyaev 1998 ). The difference is that that paper considered the 
scattering of a Planckian spectrumQ 

For a given /t, the photon frequency change, m, ranges in [(1 -/3)/(l -/3/i),(l +/3)/(l -/3/i)]. Integrating equation ( |83| ) over the 
solid angle of emergence of the scattered photon yields the spectrum that forms as a result of scattering of monochromatic radiation 
on electrons moving at a speed of (3. For a given u, the integration limits are 
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(84) 



The result is 
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(85) 



^A formula similar to equation mA) made it possible to find relativistic corrections to the amplitude of CMB distortions in the direction of clusters of galaxies. 
These corrections are of the order oijKTc/mcC^) , (V/c)^, and (V^/c) X (kT^/m^c^), where V and Vr are the cluster peculiar velocity and its component along the line 
of sight, respectively. 
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E'/E 

Fig. 17. — Spectra resulting from the single scattering of isotropic monochi'omatic radiation on an ensemble of electrons moving isotropically at a given speed, 
s calculated from formula ( p5[ ) for different values of the electron velocity (v = /3c). Compton recoil is ignored, which coiTesponds to the Hmit hv <^ ntcV^. 

where the subscript plus sign corresponds to the values m > 1, and the minus sign to m < 1. Equation (pa) was earlier derived by 



Fargionet al. (1997). The derivation of these authors differs from the one above in the order of integrations: they first considered the 



scattering of a beam of photons by a beam of electrons with a given angle between the two, and then implemented the integration 

over this angle. 

Let us now make a transition from a single electron to an ensemble of electrons with density A^e and introduce a quantity, related to 

the result of equation (p5|), that gives the probability of a scattering event calculated per unit dimensionless time t = ajNeCt, defined 

as 

1 dNiiy) 



P{i)^u,p) = 



ajc dtdv 
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Equation (Bq) describes the kernel of the integral kinetic equation for the problem of Comptonization on monoenergetic electrons. 
Let us mention its basic properties: (1) P{v ~^i>,(3) = {v/vfP{v -> v,l3), (2) j^"""Pdv = 1, and (3) j'^"""{v/v- \)Pdv = Aix- 1)/3, 
where v„,i„ = (1 -/3)/(l + /3), v,nax = (1 + /3)/(l -/?) are the minimal and maximal possible values of the photon frequency after 
scattering. Property 1 ensures that the detailed balance principle is obeyed (see also eq. [El]] for the case of Maxwellian electrons). 
In property 2, photon conservation manifests itself. Property 3 is a well-known relation (see, e.g., eq. [2.33] in the review by 
Pozdnyakov et al. 1983), which describes the rate at which the electrons transfer energy to the photons. There exists a more general 
expression than equation (pSb (which was obtained in the Thomson limit) for the Piv —^ v, f3) kernel (see eqs. [8.1 .8] and [8. 1 . 12] in 
Nagirner & Poutanen 1994, and references therein), which is valid for arbitrary photon and electron energies. 

In Figure 17, we present examples of spectra described by equations (pSb and (p6|). Their characteristic feature is the presence 
of a cusp at ly/D = 1. Note that the right wing of the line contains more photons (w 1/2 + 69/3/140) than the left one (« 1/2- 
69/3/140), the asymmetry becoming more pronounced as the electrons get more relativistic. Such spectra will build up when isotropic 
monochromatic radiation is scattered off an optically thin cloud of electrons that are movin g isotropically at the same speed. Note 
that equation (H) retains valid for ultrarelativistic (7 ^ 1) electrons (the case that interested |^argion et al . 1997 ). 

In the nonrelativistic limit (/3 <C 1), the frequency is changed upon scattering by a small amount: |m-1|/(m+1) ^ 1. This allows 
us, by expanding equation ( p5[ ) in powers of /3 and 

u— 1 
^=—7, (87) 



u+l 



to derive the formula 



P(i}^,y,p)= — 



11 
20" 



73 
140 



p^ + 0((3*)\ + 



lei 



-l^\p'^oip^) 



NARROW LINE AFTER SINGLE SCATTERING 



29 



P^ + 0((3*) 



-(3' + 0m + 



.1 + 11^2 + 0(^4) 

10 lO'^ ^ 



51 

"to 



(3^ + 0(13*) 



(1+0, 



(88) 



which approximates well the scattered photon distribution. The series given in equation (K8h ensures photon conservation to an 
accuracy of ^ Pdv = \ + 0([3''), or, putting the same expression in explicit form, ^ Pdv= 1 + 19/5'*/300+ .... 

The spectrum that forms as a result of the single scattering of isotropic monochromatic radiation on a group of electrons with a 
given isotropic distribution of velocities, fp, is described in general by the formula 



(89) 



P{v^v)= P{v^v,l3)fpl3'dl3, 



where P(£' — > z^, /?) is governed by e qua tions (pa) and (pq). In the case of a thermal plasma, one must insert the relativistic Maxwellian 
distribution function into equation (K9b. The result can then be evaluated numerically. 

In the present study, we are interested in the mildly relativistic case, ?/ = kTg/mgC^ ^ Ol- In this limit, we can derive P(t> — > u) 
analytically by making use of the approximate equation (|8^. To this end, the fp function must be expanded in terms of /?. In order 
to obtain the result with the accuracy we need, it is enough to retain only relativistic correction terms of order rj in this series, i.e., to 
substitute equation ( p^ ) of §3 for fp. 

The integral in equation (p9|) is to be performed in the range of values /? > f3,„ = \£^\ = |(m-1)/(m+1)|. Upon elementary calculations, 
we arrive at equation (EQ) (with the transition from v, v to v, v'), which was our aim. 
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